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ABSTRACT 

We obtain contemporaneous observations of the surface and corona of AB Dor, a 
^ • young single cool star, using ground-based circularly polarised spectra from the Anglo- 

&\ | Australian Telescope and X-ray lightcurves and spectra from the Chandra satellite. 

The ground-based data are used to construct surface magnetic field maps, which are 
^^O 1 extrapolated to produce detailed models of the quiescent corona. The X-ray data serve 

as a new test for the validity of these coronal models. 

We find that AB Dor's X-ray corona must be concentrated close to its surface, 
with a height, H ~ 0.3-0.4R*; this height is determined by the high coronal density 
and complex multi-polar magnetic field from the surface maps. There is also signif- 
• icant correlation between the positions of surface and coronal active longitudes as 

Oh" determined from the surface spot and magnetic field maps and the X-ray lightcurve. 

q . At this epoch (2002 December) AB Dor appears to possess one very large active longi- 

' tude region, covering almost half the star; displaying enhanced activity in the form of 

! large dark spots, strong magnetic fields and chromospheric emission. This is unusual 

C$ • as previous surface maps of AB Dor typically display more active regions that span a 

wider range of longitudes. 

Finally, the level of rotational modulation and shape of the X-ray lightcurve de- 
■ pend on the distribution of magnetic field in the obscured hemisphere (AB Dor is 

j inclined by 60°). The models that best reproduce the rotational modulation observed 

in the contemporaneous Chandra X-ray lightcurve and spectra require the magnetic 
field in the obscured hemisphere to be of the same polarity as that in the observed 
hemisphere. The Sun shows different behaviour, with the leading polarity reversed in 
the opposite hemisphere. The X-ray observations provide a unique constraint on the 
magnetic structure in the obscured hemisphere. 

Key words: stars: spots - stars: magnetic fields - stars: imaging - stars: activity - 
stars: coronae - X-rays: stars 



1 INTRODUCTION 

Magnetic activity in rapidly rotating cool stars leaves dis- 
tinctive signatures at all atmospheric levels: from stellar sur- 
faces where strong magnetic fields cause the formation of 
large dark spots, to their coronae which typically have sig- 
nificantly larger temperatures, pressures and densities than 
those seen in the solar corona. 
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donati@ast.obs-mip.fr; bhouse@head. cfa.harvard.edu; njd2@st- 
andrews.ac.uk; kw25@st-andrews.ac.uk; dupree@cfa.harvard.edu; 
acc4@st-andrews.ac.uk; Fabio.Favata@rssd.esa.int 



Early X-ray missions, ROSAT and Einstein, established 
the global properties of X-ray emitting coronae in cool stars, 
revealing that all cool stars with outer convection zones pos- 
sess hot MK coronae, with the rapid rotators displaying X- 
ray luminosities that are over 2 orders of magnitude larger 
than those observed on the Sun (Schmitt 1997). The lat- 
est X-ray and EUV missions such as EUVE, Chandra and 
XMM- Newton provide access to line diagnostics with high 
spectral resolution. This higher spectral resolution has been 
exploited to (a) probe the temperature distribution of hot, > 
MK, plasma in active stellar coronae, revealing the distribu- 
tion of coronal temperatures (e.g., Sanz Forcada, Brickhouse 
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& Dupree 2003; also see review by Dupree 2002), and to (b) 
measure coronal densities that are greater than 10 10 cm~ 3 in 
the hot MK plasma (e.g., Testa, Drake & Peres 2004, Ness et 
al. 2002, Giidel et al. 2001). These densities suggest typically 
small X-ray coronal structures (Dupree et al. 1993). 

While the physical properties of active stellar coronae 
are increasingly well determined, the detailed structure of 
stellar coronae is still a matter of debate; we have yet to 
ascertain how far the corona extends and how it is dis- 
tributed around the star. This is important as the struc- 
ture of the stellar coronal field affects how zero-age main- 
sequence (ZAMS) stars spin down (Solanki et al. 1997), how 
binary systems interact, and even how planets form (Feigel- 
son, Garmire & Pravdo 2002). Radio studies have found 
that stable magnetic fields can exist out to several stellar 
radii around active RS CVn binary systems (e.g., Mutel et 
al. 1985). There is also evidence for very extended mag- 
netic fields in the main sequence dM5.5e binary system, UV 
Cet A and B, from radio observations using the Very Large 
Baseline Array (VLB A) . Images of the strongest radio emit- 
ter appeared to have two stable radio-emitting components 
caused by stable magnetic field structures at a distance of 
between 2.2 to 4R„ (Benz, Conway & Giidel 1998). These 
sizes are consistent with the sizes of the stable prominences 
commonly detected as Ha absorption transients on rapidly 
rotating cool stars (e.g., Collier Cameron & Robinson 1989, 
Collier Cameron et al. 2002). 

While active stars are too distant to resolve directly, We 
can employ indirect techniques to probe the structure of the 
hot, MK material in their coronae using X-ray and EUV in- 
strumentation. Recent studies of binary stars have enabled 
us to gain insights into where X-ray emitting regions are lo- 
cated: e.g., periodic variations in the EUV lightcurve of the 
contact binary 44i Boo strongly indicate that the emission 
originates in a small, dense, high-latitude region on the pri- 
mary star (Brickhouse & Dupree 1998). Subsequent Chan- 
dra grating spectra confirmed this picture for 44i Boo using 
lightcurves and Doppler velocity measurements (Brickhouse, 
Dupree & Young 2001). Hussain et al. (2005; hereafter, Pa- 
per I) have employed these techniques to the analysis of AB 
Dor. 

X-ray flare analyses reveal the sizes of X-ray emitting 
coronal structures through eclipse mapping techniques; e.g., 
a flare observed in the eclipsing binary Algol (B8 V+K2 IV) 
was completely eclipsed, thus strongly constraining its loca- 
tion to 0.5 R* above the "south" pole of the K star (Schmitt 
& Favata 1999; Favata & Schmitt 1999). While flares on 
active G and K-type main sequence stars appear to origi- 
nate within 1R» above the stellar surface (e.g., Mullan et 
al. 2006), a recent Chandra study has found that young (1- 
10 Myr) stellar objects in the Orion Nebula Cluster have 
flaring loops that can extend to much greater heights, up to 
5-10 R, (Favata et al. 2005). 

AB Doradus (AB Dor, HD 36705) is an active young 
(50 Myr) K0 dwarf that has recently arrived on the main se- 
quence. It rotates very rapidly (Prot = 0.51479 day, v e sini 
= 90kms _1 ), so two rotation cycles can be covered in just 
over a day. AB Dor is also bright at both optical and X- 
ray wavelengths (my ~7, logLx/Lboi ~ 30 (Vilhu & Lin- 
sky 1987). AB Dor displays signatures of magnetic activity 
at all atmospheric levels: Doppler maps have revealed the 
presence of large dark surface spots, with a large spot cover- 



ing its pole and smaller spots at lower latitudes, and strong 
kG field covering almost all observable latitudes (e.g., Do- 
nati et al. 1999, 2003); while long-term X-ray monitoring 
indicates that AB Dor flares almost once a day (Kiirster et 
al. 1997). AB Dor's emission measure distribution (EMD) 
is typical for an active cool star, whether an RSCVn binary 
or a ZAMS star (Sanz-Forcada et al. 2002, 2003). Indeed, 
AB Dor displays characteristic magnetic activity signatures 
that appear typical of rapidly rotating cool stars in general. 

There is evidence for both compact and extended struc- 
tures in the corona of AB Dor. For example, small X-ray 
flaring loops are detected in flare decay analyses of X- 
ray lightcurves obtained using BEPPOSAX (Maggio et al. 
2000). They find that the flaring region is small (H < 0.3 R») 
and that it is uneclipsed over a full rotation cycle, indicating 
that it is located near the pole of the star. A statistical anal- 
ysis of ROSAT data spanning a period of 5 years suggests 
a small level of X-ray rotational modulation at between 5- 
13% (Kiirster et al. 1997). Brandt et al. (2001) find that 
between 60% to 80% of the UV emission should originate 
near the stellar photosphere from HST/GHRS observations 
of AB Dor. Radio emission from AB Dor shows rotational 
modulation; this modulation appears to be associated with 
magnetic structures lying over the most spotted regions, 
although the exact geometry of the radio source depends 
strongly on the specific model used (Vilhu et al. 1993, Lim 
et al. 1994). 

Conversely, there is also evidence for extended struc- 
tures in AB Dor's corona, albeit from cooler material at 
transition region temperatures. For example, transition re- 
gion lines observed in FUSE and HST datasets (C IV 1548A, 
Si IV 1393A and O VI 1032A) are observed to have broad- 
ened wings in emission extending to velocities of 270km/s. 
This broadening may be caused by optically thin plasma 
(T w 10° K) at heights out to the Keplerian co-rotation 
radius, 2.6 R» (Brandt et al. 2001, Ake et al. 2000). The 
transition region line, O VI, shows modulation with rota- 
tional phase in both flux and profile. Lower flux and nar- 
rower profiles correspond to the predicted locations of open 
magnetic field lines (Dupree et al. 2006). AB Dor also shows 
evidence of extended structures, such as slingshot promi- 
nences that extend out to over 2R«; these are detected as 
fast-moving absorption transients in Balmer lines (e.g., Ha) 
(Collier Cameron & Robinson 1989, Donati et al. 1999). It 
is not clear whether or not they are in the open or closed 
corona although Jardinc & van Ballegooijen (2005) show 
how these structures can be trapped in long thin loops that 
are in the stellar wind, lying beyond the closed, X-ray emit- 
ting corona and extending out to 5R* in rapidly rotating 
stars like AB Dor. 

We obtained contemporaneous observations of AB Dor 
at both optical and X-ray wavelengths in 2002 December 
in order to investigate the relationship between surface and 
coronal active regions. Surface spot and magnetic field maps 
have been obtained using the spectro-polarimetric optical 
data from the Anglo-Australian telescope (AAT) observa- 
tions, and the stable, non-flaring X-ray emitting corona has 
been probed using the Chandra satellite. This is the first 
time that surface magnetic field and spot maps have been 
obtained contemporaneously with X-ray observations for an 
active cool star. We aim to extrapolate our surface magnetic 
field maps to construct detailed models of the X-ray emitting 
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corona of AB Dor at this epoch. We will use these models to 
predict the amount of X-ray rotational modulation that we 
expect to see in this star assuming an isothermal corona and 
compare our predictions to the Chandra X-ray observations. 



2 MULTI- WAVELENGTH CAMPAIGN OF AB 
DOR IN 2002 DECEMBER 

The active KOV star, AB Dor, was observed contemporane- 
ously using the low energy transmission grating (LETG) and 
the high resolution camera spectroscopic array (HRC-S) on 
the Chandra X-ray satellite, and the SemelPol polarimeter 
on the UCL echelle spectrograph (UCLES) on the Anglo- 
Australian Telescope (AAT). The X-ray observations lasted 
almost two full rotation cycles (commencing at 21:25 UT on 
2002 December 10 and ending at 21:53 UT on 2002 Decem- 
ber 11). No large flares were observed over this period; by 
covering almost two rotation cycles we could study X-ray 
rotational modulation caused by the quiescent non-flaring 
corona. The ground-based AAT observations spanned five 
consecutive nights (2002 December 11-December 15), over- 
lapping with the X-ray observations for eight hours. Sur- 
face magnetic fields and spots on AB Dor remain stable on 
timescales of over 1 week (e.g., Donati et al. 1997). Hence 
no significant changes will have occurred in the surface spots 
and magnetic field distributions over the course of the Chan- 
dra X-ray observations. 

The high resolution spectro-polarimetry obtained from 
the SemelPol/UCLES/AAT configuration is inverted to re- 
construct surface spot and magnetic field maps of AB Dor. 
Doppler imaging techniques were used to invert the AAT 
spectra and reconstruct detailed surface spot and magnetic 
field maps of AB Dor (see Donati et al. 2003, 1999). Full 
phase coverage was obtained over the observing run, result- 
ing in surface maps with a latitude resolution of approxi- 
mately 3°. 

2.1 Chandra X-ray observations: main findings 

The Chandra X-ray observations are described in detail in 
Paper 1. Briefly, AB Dor was observed over one continuous 
pointing spanning 88.1 ksec (Observation ID 3762), which 
corresponds to 1.98 P ro t- Paper 1 reports detailed evidence 
of rotational modulation observed in both the X-ray 
lightcurves and spectra; the main findings from Paper 1 are 
summarised in Figure 1 and listed below: 

• There is significant rotational modulation in the X-ray 
lightcurve, which has a flat level of emission superimposed 
with three peaks that cause 12% rotational modulation. 

• The O VIII 18A line profile shows a rotationally mod- 
ulated velocity pattern that repeats in two consecutive ro- 
tation cycles. The best-fit sine curve has a semi-amplitude 
of 30 ± 10 km s" 1 . 

• The Chandra HETG Fe XVII 15A line profiles indi- 
cate that the X-ray emitting corona does not extend beyond 
0.75 R» above the surface, assuming that Doppler broaden- 
ing is the cause of any excess broadening in the spectral line 
profile. We note also that FUSE observations of the coro- 
nal forbidden Fe XVIII and Fe XIX lines at 974 and 1118 A 
suggest an even smaller upper limit of approximately 0.5 R, . 



2.2 Ground-based AAT observations 

The SemelPol/UCLES/AAT configuration, used to obtain 
the ground-based spectro-polarimetric data, has been de- 
scribed in detail in previous papers (Donati et al. 2003, Do- 
nati et al. 1999, Donati & Cameron 1999). Briefly, SemelPol 
is a visitor instrument, mounted at the Cassegrain focus of 
the 3.9-m AAT and coupled with the UCLES spectrograph. 
The CCD chip records two beams of opposite polarisation 
spanning a wavelength range of over 2400A, from 4374.7A 
(order 129) to 6816. 5A (order 84). The observations of AB 
Dor were obtained between 2002 December 11-15 (see Table 
1 for details). 

The cross-correlation technique of least squares decon- 
volution (LSD) is used to sum up the signal from over 2700 
photospheric lines, hence exploiting the full wavelength cov- 
erage achieved with the AAT/SEMEL instrument configura- 
tion (Donati et al. 1999). The reduced spectra had a typical 
peak S:N level of approximately 146 to 200, the typical mul- 
tiplex gain in S:N was a factor of 40-41, and the output LSD 
S:N ranges between 5895 and 8512A (also see Table 1). 

2.3 Surface magnetic activity maps of AB Dor 

Doppler imaging (DI) techniques enable us to produce de- 
tailed surface spot and magnetic field maps of magnetic ac- 
tivity in rapidly rotating cool stars by inverting a time-series 
of high resolution spectra. The rotational broadening in the 
line profiles of active rapidly rotating cool stars separates 
the distortions caused by surface spots and magnetic fields 
in velocity-space (Vogt, Penrod & Hatzes 1987, Semel 1989, 
Collier Cameron & Unruh 1994, also see review by Hussain 
2004). The DI code, DoTS, was used to reconstruct the spot 
map of AB Dor (Figuref2| by inverting a time-series of inten- 
sity line profiles (Collier Cameron 1997). The technique of 
Zeeman Doppler imaging (ZDI) enables the reconstruction 
of surface magnetic field maps by inverting a time-series of 
circularly polarised spectra from rapidly rotating stars. The 
code described by Donati & Brown (1997) was used to pro- 
duce the surface magnetic field maps shown in Figure[2] 

AB Dor's spot distribution has been mapped for over 
a decade. Its spot patterns are typical of rapidly rotating 
active stars in general, with a large polar spot covering the 
polar region and extending to below 70° latitude. This polar 
spot appears to be almost unchanged for at least seven years. 
ZDI is a more recent technique, hence until recently AB 
Dor was the only active rapidly rotating star for which we 
had detailed magnetic field maps (latitude resolution ~ 3°). 
The magnetic field maps in 2002 December are also very 
typical for AB Dor, with strong field covering almost all 
observable latitudes. Strong azimuthal field is found, with 
a strong unidirectional component at high latitudes; this 
is very similar to previous magnetic field maps of the star 
(Donati & Collier Cameron 1997; Donati et al. 1999; Donati 
et al. 2003). 

2.4 Comparison of active longitudes at the 
surface and corona 

The data acquired here enable us to compare the longitude 
positions of active regions at all atmospheric levels: from 
the surface (as probed by the surface spot and magnetic field 



4 G. A. J. Hussain, et al. 




0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

Rotation phase Rotation phase 



Figure 1. (a) The X-ray lightcurve shows significant rotational modulation at the ±12% level. The lightcurve shows three components 
of emission: (i) a flat unmodulated level of emission (solid horizontal line), (ii) three peaks (emphasised here by outlining the peaks with 
thick solid lines) that appear to repeat in the two rotation cycles and thus are rotationally modulated, and (iii) intrinsic flaring variability 
on short timescales. (b) The centroid of the OVIII 18A line profile traces (±30 kms -1 ) velocity shifts that appear to display rotational 
modulation over two consecutive rotation cycles. See Paper 1 for more details. 



Table 1. Observing log for 2002 December 11-15. The UT dates and 
times for each exposure are shown in the first three columns; n eX p is 
the number of intensity spectra/polarised spectra taken over the times 
listed; and t eX p is the the exposure time. The first row is the Chandra 
observation, and the remaining rows are the the AAT observations. 



Date 


MJD 


UT (h:m:s) 


ricxp 


tcxp 


2002 Dec. 10-11 


52618-9 


21:25-21:53 




88100 


2002 Dec. 11 


52619.0 


10:30:40/11:22:09 


12/3 


200 






11:50:53/12:41:24 


12/3 


200 






13:13:19/14:22:16 


16/4 


200 






14:49:58/15:59:01 


16/4 


200 






16:28:46/17:37:41 


16/4 


200 






18:07:54/18:27:52 


4/1 


200 


2002 Dec 12 


52620.0 


09:57:15/11:06:16 


16/4 


200 






11:39:42/12:48:41 


16/4 


200 






13:25:19/14:34:26 


16/4 


200 






15:01:52/16:18:07 


16/4 


200 






16:45:34/17:17:52 


8/2 


200 






17:45:55/18:18:12 


8/2 


200 


2002 Dec 13 


52621.0 


09:45:13/10:35:51 


12/3 


200 






11:02:54/11:53:34 


12/3 


200 






12:22:15/13:12:52 


12/3 


200 






13:39:23/14:30:02 


12/3 


200 






14:57:26/15:48:06 


12/3 


200 






16:22:31/17:13:05 


12/3 


200 






17:40:35/18:12:44 


8/2 


200 


2002 Dec 14 


52622.0 


09:37:02/10:34:30 


12/3 


200 






10:54:35/11:45:11 


12/3 


200 






15:49:18/17:35:04 


24/6 


200 






18:03:31/18:17:20 


4/1 


200 


2002 Dec 15 


52623.0 


11:08:20/11:40:36 


12/3 


200 



maps) to the corona (as probed by the X-ray lightcurve). As 
discussed in Paper 1, the X-ray lightcurve can be described 
as a flat unmodulated emission level superimposed with the 
enhanced X-ray emission from two or three active regions; 
these active regions cause peaks in the lightcurve near phases 
0.2, 0.5, and 0.9. The peaks also seem to be associated with 
short-lived energetic flaring regions. The lower envelope of 
emission repeats fairly well in both cycles (particularly the 



peaks near phases 0.9 and 0.5 ), suggesting the presence of 
coronal X-ray emitting structures that remain stable for the 
period of observation. Figure^ shows the X-ray lightcurve: 
this has been normalised using the peak counts in the stable, 
rotationally modulating part of the X-ray lightcurve. We 
attribute the emission above 1.0 in this plot to emission 
from short-lived flaring regions. 

Figures[3jb-d) have been computed using the surface 
spot and magnetic field maps in Figure]!?] Spot coverage 
(f s ) ranges from 1.0 (representing spots covering the en- 
tire projected stellar disk) to 0.0 (representing unspotted 
photosphere). Figure[3jd) was computed by projecting the 
surface spot map onto the inclined stellar disk at each rota- 
tion phase covering a full rotation cycle. The spot coverage 
at each rotation phase is then computed by integrating the 
spot coverage over the entire inclined stellar disk. A similar 
process was carried out with the surface magnetic field maps 
to compute B 2 and -Bproj as a function of rotation phase. 
Figure|3p shows the square of the line-of-sight component 
of the magnetic field (B pro j = |Bsinf)| 2 ) as a function of 
rotation phase. Figures[3}5-d show the presence of a large 
active region at the stellar surface that peaks near rotation 
phase 1.3, very near the centre of the two largest peaks in 
the X-ray lightcurve (Fig|3K). 

The significance of this correspondence can be quan- 
tified using the non-parametric statistical Kendall r test. 
This test measures the strength of dependence between two 
variables, quantifying the degree of correlation or lack of 
correlation. We investigate whether or not there is any cor- 
relation between the positions of surface spots and X-ray 
emission in AB Dor. As required by the Kendall test, both 
datasets need to have the same binning. We adopt the same 
time-bins used in the X-ray lightcurve when computing the 
f s and \B 2 \ vs. phase plots; the observation phases between 
0.10 to 0.23 are omitted in the subsequent calculation as 
there are phase gaps in this region in the surface maps and 
therefore spot and magnetic field detections in these areas 
are less reliable. Similarly the flaring regions (above 1.0 in 
the top plot of Figure[3]) are omitted in these calculations 
as we are only interested in the quiescent coronal emission. 
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Figure 2. Surface maps of AB Dor derived using SemelPol/UCLES/AAT spectra. These maps are rectangular projections of the stellar 
surface, latitudes below -30° are not observed as the star is inclined at 60° . Observation phases are denoted as tick marks along the top 
plot. Top: Spot map of AB Dor from 2002 December; black denotes full spot coverage and white denotes unspotted photosphere. Middle 
and Bottom: Radial and azimuthal magnetic field maps of the stellar surface; blue and red denote +/- 1300 G. These surface spot and 
magnetic field maps are characteristic of AB Dor; with strong field covering all observable latitudes and the presence of a large polar cap 
in the spot map. The fits to these data will be presented by Donati et al. (in prep.). 
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Table 2. r test of magnetically active regions 
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Figure 3. Comparison of active longitudes in coronal (X-ray) 
and surface (magnetic field and spots) active regions, (a) Chan- 
dra X-ray lightcurve — asterisks and filled symbols represent the 
first and second rotation cycle respectively. The lightcurve has 
been normalised so that the lower envelope of the X-ray emis- 
sion, representing the component that is rotationally modulated 
lies below 1.0. X-ray emission above 1.0 is caused by short-term, 
most likely flaring, variability and is excluded from our calcula- 
tions for the Kendall r-test. Figures (b)-(d) have been computed 
using the surface magnetic field and spot maps shown in Fig- 
ure 2; (b) Bproj is the total magnetic field projected along the 
line-of-sight summed over each phase; (c) the total magnetic field 
squared, B 2 from the Zeeman Doppler maps shown in Figure 2; 
and (d) the spot coverage (/ s ) computed from the spot map in 
Figure 2 as a function of rotation phase. 



The r test is used to rank-order the f s and X-ray emission, 
and the \B 2 \ and X-ray emission; it enables us to measure 
the degree of correlation, the statistical probability of this 
correlation being due to chance (Press et al. 1988). We find 
highly significant correlation between all three quantities in- 
vestigated: i.e. X-ray flux, spot coverage and \Bg\. 



3 BUILDING A CORONAL MODEL OF AB 
DOR 

Jardine et al. (2002) extrapolate surface magnetic field maps 
of AB Dor and use the resulting three-dimensional magnetic 



field models to predict coronal properties of rapidly rotat- 
ing active stars (e.g., X-ray emission levels and rotational 
modulation, locations of open field; Jardine et al, 2006). 
The contemporaneous AAT and Chandra datasets enable 
us to test our X-ray modelling capabilities in detail for the 
first time. The X-ray corona is modelled using the following 
inputs: (a) the three-dimensional magnetic field model ex- 
trapolated from the surface magnetic field map (Figure 3), 
(b) an isothermal corona that has a temperature of approx- 
imately 10 MK (consistent with emission measure distribu- 
tions of AB Dor), and (c) assuming hydrostatic equilibrium. 
The extent of the closed corona depends on the strength and 
complexity of the surface magnetic field maps, the emission 
measure level and the electron density of hte X-ray corona 
of the star (as determined from X-ray observations of AB 
Dor). This dataset represents a unique opportunity to use 
measurements from X-ray data obtained contemporaneously 
with magnetic field maps to test our models and help in con- 
structing more realistic models. 

3.1 Coronal fluxes and densities in 2002 
December 

The Chandra data were re-reduced using the latest dedicated 
CIAO pipeline processing routines (CIAO v. 3. 3. 0.1). The 
most prominent He-like triplet, the OVII triplet, is used 
to to measure AB Dor's coronal density in 2002 December. 
The resonance, intercombination and forbidden lines are all 
detected clearly in these X-ray spectra. Measurements based 
on previous observations (Testa, Drake & Peres 2004, Sanz- 
Forcada, Maggio & Micela 2003) can be compared with our 
measurements to estimate the stability of the densities in 
the X-ray emitting regions of the stellar corona. 

In order to measure line fluxes for this triplet we first 
had to evaluate the continuum level accurately. We estimate 
the continuum level by generating fake or model continuum 
and spectra using the EMD derived by Sanz-Forcada, Mag- 
gio & Micela (2003) from Chandra/HETG /ACIS-S spectra 
of AB Dor obtained in 1999 October (Observation ID 16). 
The spectra are folded through the LETG instrument re- 
sponse function, modelling for combined orders up to order 
5. We find that the fluxes from these model spectra are a 
good fit to the observed spectra with very little scaling re- 
quired suggesting that the emission measure distributions 
are very similar at both epochs (1999 October and 2002 De- 
cember; Sanz-Forcada, Maggio & Micela 2003). These model 
spectra are used to determine the best line-free regions, and 
the model continuum is scaled to fit these areas. As these 
regions are relatively small, there is still considerable uncer- 
tainty in the placing of the continuum level and this con- 
tributes the largest uncertainty to the subsequent flux mea- 
surements. Once the continuum level is determined, each 
line profile is fitted using a least-squares Gaussian fitting 
procedure; the exact shape of the profile fit does not make 
much difference to the final measurements. We also fit the 
Moffat function, i" = (1 + (A/A c ) 2 ) _/3 , with f3 = 2.5, as this 
shape is a closer representation of the instrumental line pro- 
file (Chandra Proposers' Observatory Guide Section 9.3.3 
v.8), however the resulting fluxes are consistent with the 
Gaussian fitting method. The measured fluxes are reported 
in Table 3 and the fits are shown in Figure[4] 

The ratios: f/i = 2.1 ± 0.2 and (/ +i)/r = 0.86 ± 0.1 
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Table 4. Flux measurements for the He-like OVII triplet. 



Wavelength (A) 21.6 (r) 21.8 (i) 22.1 (f) 

Fluxes (10 _5 )erg/cm 2 /s 26.8 7.4 15.6 

Errors (10 _5 )erg/cm 2 /s 1.9 1.3 1.6 



100 - 



/I i 




21.4 21.6 21.0 22.0 22.2 

Wavelength 

Figure 4. The He-like density-sensitive OVII triplet at 2lA. 



correspond to an electron density, logn e [cm -3 ] ~ 10.3to'3- 
The (/ + i)/r ratio gives the electron temperature, and is 
found to be 2.2 MK (using the Astrophysical Plasma Emis- 
sion Database, APED; Smith et al. 2001). The errors quoted 
for the density take the largest source of uncertainty into ac- 
count; i.e. the placement of the continuum level (Brickhouse 
2002). This uncertainty dominates the measurement errors 
associated with line profile fitting. For instance, if the con- 
tinuum is actually lower than the value we adopt (e.g., scaled 
by 80% of the current value), this would increase the mea- 
sured n e by a factor of 2, while still fitting the continuum 
level in the line-free regions of the X-ray dataset. 

Measurements of densities at higher coronal tempera- 
tures are harder to determine with our LETG dataset as 
either the diagnostics lines are at short wavelengths where 
the spectral resolution is poor or else they are at very high 
wavelengths (91-128A) where the background level is high 
and the lines are relatively weak. 

AB Dor's coronal n e has previously been estimated us- 
ing the OVII He-like density sensitive ratio from the data 
obtained using MEG/ACIS-S and XMM-Aeu>fon/RGS ob- 
servations of the star. A full table of previous measurements 
of n e are shown in Table [3~T1 These measurements show that 
our measured value of n e is consistent with previously mea- 
sured values, especially given the error introduced by the 
difficulties of measuring the continuum level accurately. The 
NelX and MgXI diagnostics, formed at higher tempera- 
tures, show higher densities. Using Chandra/HETG spec- 
tra, Sanz-Forcada et al. (2003) report logn e [cm~ 3 ] ~ 11 
using NelX (T = 3.98 MK) and logn e [cm~ 3 ] > 12; Ness et 
al. (2004) measure log n e [cm~ 3 ] = 11.4 ± 0.16 using NelX 
(T = 4MK); and Testa et al. (2004) measure log7i e [cm~ 3 ] < 
12.25 using MgXI (T = 3.98MK). 



3.2 Constructing models: reconciling field, EMD 
and n e measurements 

We build a simple isothermal coronal X-ray model by ex- 
trapolating the surface magnetic field maps in Figure]!?] All 
of the closed field lines are assumed to be filled with isother- 
mal X-ray emitting gas in hydrostatic equilibrium. We take 
the temperature of the corona to be T ~ 10 MK in accor- 
dance with the stable peak in AB Dor's EMD (e.g., Sanz- 
Forcada et al. 2003). Before we compute the X-ray proper- 
ties of the corona at this epoch we have to determine two 
parameters: the radius of the source surface - this deter- 
mines the extent of the closed, X-ray emitting corona and 
is the point beyond which all the field becomes radial, and 
the pressure at the base of the X-ray emitting loops. The 
pressure at the footpoint of closed loops, po, scales with 
the magnetic field strength, B, by a factor K: po = KBq. 
This scaling factor, K, and the extent of the closed, emit- 
ting corona depend on the following parameters: magnetic 
field strengths from the surface field maps, coronal den- 
sity, temperature and emission measure, EM. The X-ray 
values we are trying to fit here are: coronal temperature 
(T C =10MK), density (log n e [cm -3 ] = 10.3) and emission 
measure logEM[cm -3 ] = 52. 

The coronal X-ray emissivity is assumed to be propor- 
tional to the density squared and we use a three-dimensional 
Monte Carlo radiation transfer code to produce optically 
thin X-ray images (assuming no scattering or absorption of 
the coronal X-rays). Figure[S] shows contours of goodness-of- 
fit to the observed X-ray parameters (i.e. T c , n e and EM) of 
our extrapolated magnetic field model as a function of the 
radius of the source surface, Rss, and pressure scaling factor, 
K. As the coronal models are in hydrostatic equilibrium, the 
model densities are actually an emission-measure weighted 
average over the volume of the corona. Figure[S] shows the 
dependence of the extent of the X-ray emitting corona on 
the coronal density; unsurprisingly low densities are consis- 
tent with more extended coronae. The extent of the corona 
can also be increased by increasing magnetic field strengths 
and/or increasing the EM value or temperature. However, 
given AB Dor's surface magnetic field map obtained using 
ZDI, and its coronal EM, n e and temperature, our results 
suggest an X-ray emitting corona that is very compact with 
a source surface at radius, R SB = 1.4 R*. 



4 TESTING THE CORONAL MODELS 

One of the main objectives of this paper is to learn whether 
or not these magnetic field maps and X-ray observations can 
be used to build more realistic detailed X-ray models. As AB 
Dor is inclined at 60°, we cannot observe the magnetic field 
in the obscured hemisphere. We model this "missing" field 
by reflecting the magnetic field in the observed hemisphere 
into the obscured hemisphere. Figure[6] shows the resulting 
X-ray models for two cases: when the observed field is sim- 
ply reflected into the lower, obscured hemisphere; and when 
the observed field is reflected and its polarity is reversed 
in the lower hemisphere. As shown in Figure(6] the differ- 
ent types of model do not affect the X-ray emission from 
the observed upper hemisphere; however, the reversed po- 
larity model causes more X-ray emission to originate near 
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Table 3. Density measurements for the He-like OVII triplet. 



Density Temperature Source Date Reference 

log n e [cm~ 3 ] [MK] 



10.78 ± 0.14 


1 6 - 9 
1 - D -0.5 


MEG 


1999 Oct 


Sanz-Forcada et al. 


2003 


10.2 ±0.74 


2 


MEG 


1999 Oct 


Ness et al. 2004 




10.75°J5 12 


2 


MEG 


1999 Oct 


Testa et al. 2004 






2 


RGS 


2000 May 


Giidel et al. 2001 




10.5 ±0.10 


2 


RGS 


2000 May 


Ness et al. 2004 




10.77 ±0.08 


1.7 ± 0.3 


RGS 


2000 Jun 


Sanz-Forcada et al. 


2003 


10.82 ± 0.06 


2.1± 0.3 


RGS 


2001 Jan 


Sanz-Forcada et al. 


2003 


10 3 - 4 


2.2 


LETG 


2002 Dec 


this paper 






10.8 



Figure 5. Left: Contours of goodness-of-fit that are used to select the best values of the radius of the source surface, R BB and, K, the 
scaling factor determining the pressure at the base of X-ray emitting loops. Right: Higher densities necessitate a more compact corona. 
This plot shows the relationship between the extent of the X-ray emitting corona and the coronal density, assuming the magnetic field 
map shown in Figure 1 and an isothermal corona with T=10 MK and a log (EM)[cm~ 3 ] = 52. Our measured density, logn e [cm — 3 ] = 10.3, 
suggests a very compact corona with a source surface at a radius, R BB ~ 1.4R, 



the equator of the star as more field lines are introduced that 
connect the upper and lower hemispheres if the polarities are 
reversed filling in some areas. 

Two sets of models are generated following the above 
procedure. In group A: the field in the observed hemisphere 
is reflected into the lower hemisphere, and this lower hemi- 
sphere field is then offset relative to the field in the ob- 
served hemisphere by different amounts in longitude-space; 
X-ray lightcurves are computed for each of these models 
and the rotational modulation evaluated. In group B: the 
field is reflected, the polarity is reversed and this modi- 
fied field is then displaced relative to the field in the up- 
per hemisphere by different amounts in longitude; as with 
group A, X-ray lightcurves are computed for each model and 
rotational modulation measured. The observed level of rota- 
tional modulation in the lightcurve determined in Paper 1 is 
±12%. We verified that this level of rotational modulation 
is reliable by also constructing a X-ray lightcurve using just 
the flux in the strongest line profiles, hence reducing any 
contribution from the continuum flux in the dataset. The 
resulting lightcurve also shows a consistent level of ±12% 
rotational modulation superimposed with stochastic short- 
term variability. 



4.1 X-ray lightcurves 

Paper I uses a sine-curve fitting technique to quantify the 
rotational modulation in the Chandra X-ray lightcurve and 
measure ± 12% rotational modulation from the mean level 
of the lightcurve. This is consistent with previous estimates 
of X-ray rotational modulation for AB Dor, with estimates 
varying between ±5-13% (Kiirster et al. 1997; Brandt et al. 
2001). 

Lightcurves generated from our two sets of models are 
shown in Figure[H] These plots confirm that the field in the 
obscured hemisphere would have a significant effect on the 
observed X-ray rotational modulation. The models in group 
A show rotational modulation of approximately ±8%, while 
models in group B show much larger levels of X-ray rota- 
tional modulation (±21%). 

We find that the detailed shape of the lightcurves can 
vary considerably in models with different phase offsets be- 
tween the two hemispheres. However, the overall amount of 
rotational modulation in the models in each group does not 
vary significantly. These findings suggest that the detailed 
shape of an observed X-ray lightcurve would be hard to fit 
as it is strongly dependent on field that we cannot detect di- 
rectly in ZDI surface magnetic field maps; we therefore con- 
centrate on understanding the family of models for which we 
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Figure 6. Snapshots of the spot map and the X-ray models at rotation phases 0.3 (left column) and 0.8 (right column). Top: "Snapshots" 
of the surface spot distribution on AB Dor in 2002 December: grey and black represent unspotted and spotted regions of the photosphere 
respectively. The top left and top right images show the star at the maximum and minimum spottedness phases respectively. Middle 
row: X-ray images generated from a magnetic field model in which the observed field is directly reflected into the obscured hemisphere 
(as the star has an inclination angle of 60°). This is called model 'a'. Bottom left and right: In this X-ray model, the observed field is 
reflected into the obscured hemisphere and its polarity is reversed; this is called model 'b'. The X-ray emission pattern remains similar in 
the top hemisphere in both the 'a' and 'b' X-ray models, the main difference in model 'b' is that there are more field lines connecting the 
top and bottom hemispheres near the equator, thus increasing the X-ray emission contribution from this region. The X-ray lightcurves 
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Figure 7. Left: X-ray lightcurve for model 'a' in Figure|6] Right: X-ray lightcurve for model 'b' in Figure|6] Model 'b' has a much higher 
level of X-ray rotational modulation than observed for AB Dor at this, or any other, epoch; therefore it is likely that the field in the 
obscured hemisphere is of the same polarity as that in the opposite hemisphere. 



can reconstruct the same level of X-ray rotational modula- 
tion as that observed in the Chandra dataset. The group A 
set of models provides the best fit to the level of rotational 
modulation observed in AB Dor, not only at this epoch but 
also at the previous epochs for which there have been long 
X-ray observations of the star (Kiirster et al. 1997). This 
suggests that the field in the obscured hemisphere is of the 
same polarity as that in the magnetic field map, unlike the 
solar case, where the leading polarities are reversed in op- 
posite hemispheres. 

4.2 X-ray line profiles 

Using the coronal X-ray models shown in the previous sec- 
tion, we produce simulated X-ray line profiles assuming the 
corona is optically thin and co-rotating with the star. In 
forming the X-ray images, we now take into account the 
Doppler shifts of the photons emitted from different parts 
of the corona. We assume the X-ray line emission is at a sin- 
gle rest-frame wavelength (infinitely sharp line emission) and 
by binning the emitted photons according to their Doppler 
shifts we form X-ray channel maps. Line profiles are then 
formed by integrating the total intensity in the individual 
channel maps. 

Forty line profiles are produced, evenly sampling a full 
rotation cycle. We phase-bin the line profiles to produce 
quarter-rotation-phase binned line profiles for three different 
X-ray models: (a) with a source surface at 1.3 R* and the 
magnetic field reflected into the obscured hemisphere, (b) 
with a source surface at 1.3 R» with the reflected magnetic 
field in the obscured hemisphere possessing the reversed po- 
larity and (c) with the magnetic field reflected into the ob- 
scured hemisphere and a source surface at 3 R* . This last 
model has the same temperature and density as the other 
two models but with a higher EM (log EM [cm -3 ] > 53). 
The line profiles are then convolved with an instrumental 
profile corresponding to the 18.9 A line in the LETG. We 
measure the velocity shifts for these models using the same 
method as those found in Paper 1 (also shown in Figure 1), 
i.e. by fitting Gaussian profiles in the data analysis package, 
Sherpa. 

The velocity shifts found in each model depend some- 



what on the exact timing used, choosing a slightly different 
starting point can change the maximum size of the veloc- 
ity shifts observed. However, the model (a), which was con- 
structed by reflecting the magnetic field into the obscured 
hemisphere produces the largest velocity shifts (Figure |5J). 
Models (b) and (c), which possess reversed polarities and 
an extended source surface radius respectively, show signif- 
icantly smaller velocity shifts than those measured in the 
Chandra spectra. This is probably because they posses more 
uniformly distributed X-ray emission that effectively fills in 
the centres of the line profiles and reduces the measured 
velocity shifts. 

As expected, increasing the spectral resolution does 
not affect the centroid measurements significantly provided 
there is sufficient signal: increasing spectral resolution by 
a factor of four alters the positions of the line centroids 
by ±1.5 km/s. Naturally, improving the time-resolution, and 
producing more phase bins over the rotation cycle would in- 
crease the sizes of the velocity shifts observed. If the line pro- 
files were of higher S:N, then one would also observe larger 
velocity shifts in the peaks of the line profiles. With optical 
Doppler imaging the line profiles are inverted to obtain sur- 
face maps; these X-ray line profiles have significantly lower 
spectral and time resolution and so a similar inversion would 
not be as effective here. These results suggest that forward- 
modelling: creating X-ray models from surface field extrap- 
olations and testing these using the rotational modulation 
observed in both X-ray lightcurves and spectra is a more ef- 
fective approach in the case of lower resolution spectra such 
as these. 

We use our models to test whether or not we can dis- 
tinguish between X-ray coronae with source surfaces at dis- 
tinctly different heights using X-ray spectral line profiles. In 
Paper 1, we use the Fe XVII 15.0 A line profile from Chandra 
High Energy Grating (HEG) archive observations of AB Dor 
to measure the line profile width of the X-ray emission model 
and estimate that the bulk of the emitting corona does not 
extend beyond 1.3 R*. We simulate this using our models 
by producing line profiles at fourty evenly spaced phases, 
integrating these over the full rotation cycle, and convolv- 
ing the integrated line profiles with the instrumental profile 
for the FeXVII 15 A line in the Chandra high energy grat- 
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Figure 8. Left: X-ray lightcurves for the sets of models in which the observed field is reflected into the obscured hemisphere. Introducing 
different offsets between the observed field and the reflected field changes the structure of the lightcurve but the level of rotational 
modulation remains the same (±8%). Right: X-ray lightcurves for the sets of models in which the observed field is reflected into the 
obscured hemisphere and its polarity switched. As with the first set of models, different offsets between the two hemispheres affects the 
structure of the lightcurve but not the level of rotational modulation observed. 
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Figure 9. The error bars and crosses show the velocity shifts mea- 
sured in the line ccntroids of the observed LETG OVIII 18. 9A pro- 
files. The asterisks show the velocity shifts measured for our sim- 
ulated profiles for the model that best reproduces the lightcurve 
modulation level observed. 



ing (HEG) respectively. However, as shown in Figure[l0] by 
convolving the integrated spectral line with the instrumen- 
tal profile of the HEG FeXVII line, the current instrumental 
resolution is insufficient to distinguish between X-ray coro- 
nae that extend to 1.3 R, and 3 R» . If the spectral resolution 
is improved by a factor of two, (i.e. if the instrumental profile 
has a FWHM~ 80km/s) then the effect of extended X-ray 
emission would be visible as excess emission in the wings of 
the 3R» model, although very high S:N would be required 
to distinguish between different coronal extents as the flux 
in the broadened wings is still only about 10% of the peak 
flux in the line profile. The main evidence for AB Dor's com- 
pact X-ray corona (H ~ 1.4 R*) therefore comes from the 
mixed polarities measured on its surface combined with the 
very high electron densities measured in the Chandra X-ray 
dataset. 
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Figure 10. Synthetic spectral line profiles produced by integrat- 
ing the simulated X-ray line profiles from coronal models 'a' and 
'c' over a full rotation cycle and convolving with the HEG in- 
strumental profile corresponding to the HEG Fe XVII 15A line, 
(FWHM= 0.009 A). The flux has been normalised by the max- 
imum flux to enable easy comparison of both line profiles. The 
dashed line is for the compact X-ray model 'a' (R BB = 1.3 R*) 
and the solid line is the extended model 'c' (i? sa = 3R»). The in- 
strumental resolution is insufficient to distinguish between X-ray 
coronae with different coronal extents reliably. 



5 THE EXTENT AND FILLING FACTOR OF 
THE X-RAY CORONAL MODEL 

Measurements of pressure in AB Dor's corona are extremely 
high, over two orders of magnitude larger than those found 
on the quiet Sun. Data from HST, EUVE, Chandra, and 
XMM- Newton indicate that electron densities are similarly 
elevated in AB Dor, and other similarly magnetically ac- 
tive rapidly rotating cool stars. n e ~ 2 — 3 x 10 12 cm -3 , at 
transition region temperatures (about 3 x 10 4 K) (Brandt 
et al. 2001). At higher temperatures, around 3 x 10 6 K, 
recent XMM-Newton results indicate electron densities of 
n e » 3 x 10 10 cm" 3 (Giidel et al. 2001, Sanz Forcada et 



12 G. A. J. Hussain, et al. 



-0.9F 




9.8 10.0 10.2 10.4 10.6 10.8 
log( density) [cm -3 ] 

Figure 11. The dependence of the X-ray corona filling factor on 
the coronal density for a coronal model with a log EM[cm — 3 ] = 52 
and a T ~ 10 MK. The filling factor is defined in the text. 

al. 2003, Testa et al. 2004). At diagnostics formed at even 
higher temperatures: between 6-10 MK there is evidence for 
extremely high densities (10 12 < n e < 10 13 cm -3 ). Sanz- 
Forcada et al. (2003) and Testa et al. (2004) argue that 
these very high densities indicate the presence of very com- 
pact structures, with heights of under 17?*. 

The volume filling factor effectively measures the frac- 
tion of the coronal volume that is emitting in our X-ray 
model. The volume filling factors for our coronal models de- 
pend on the coronal density (Jardine et al. 2002), and are 
defined as follows: 

J nldV 

As shown in Figure lTTI / ~ 0.11 for our model that has a 
10 MK corona and a logn e [cm -3 ] = 10.3. For higher den- 
sities, logn e [cm -3 ] > 11, the filling factor would decrease 
as the extent of the corona is smaller; for these models 
this would mean the X-ray lightcurves would show more 
rotational modulation than we observed with the Chandra 
dataset. Our models cannot reconstruct structures that are 
small enough to reproduce these very high densities as they 
are extrapolated from the surface maps. The sizes of the 
smallest coronal structures that can be modelled using this 
technique is constrained by the spatial resolution limit of 
the surface maps (~ 3°). 

5.1 The effect of changing the coronal 
temperature on the sizes of X-ray emitting 
structures 

The stable peaks near 10 MK temperatures in EMDs of cool 
stars may be caused by the presence of small, dense flar- 
ing structures that are always present on active cool stars 
(Cargill & Klimchuk 2006). If this is the case, we may have 
overestimated the temperature of AB Dor's stable, quies- 
cent, non-flaring corona. We investigate here the effect of 
lowering the coronal temperature on the coronal extent. We 
have generated a series of models for a corona that is at 

2.2 MK, (i.e. the temperature determined from the density- 
sensitive line diagnostics at OVII 21 A). We use these to 
determine how changing the coronal temperature affects 



the extent of the X-ray emitting corona, _R SS . The model 
that best fits this 2.2 MK corona, and the measured den- 
sity, logn e [cm -3 ] ~ 10.3, will have a source surface ra- 
dius of 1.3 R*. It was not possible to find a solution for 
an EM of logEM[cm -3 ] = 52 and it had to be lowered 
to logEM[cm -3 ] = 51.8. This is consistent with what is 
observed, as the EMDs of AB Dor obtained from EUV and 
X-ray spectra indicate that the EM of cooler 2.2 MK plasma 
should indeed be lower than that of the 10 MK plasma. The 
resulting X-ray model indicates that changing the tempera- 
ture by this amount does not significantly affect the extent of 
the X-ray corona as the high electron densities limit the size 
of the X-ray corona (Figure [12)) . The rotational modulation 
in the lightcurves for this more compact model increases to 
±15% in this lower temperature 2.2 MK model (compared 
to ±8% in the 10 MK model). The filling factor increases 
slightly in this lower temperature model to / ~ 0.15. 



6 PROMINENCES AND PLAGES IN 
DECEMBER 2002 

The AB Dor Ha line is subject to variability caused pri- 
marily by bright surface plage regions and rapidly moving 
dark absorption features. The latter are attributed to cool 
(8,000 to 10,000 K) clouds of neutral hydrogen in enforced 
co-rotation with the stellar surface (Collier Cameron and 
Robinson 1989a). By measuring the rate at which the ab- 
sorption features drift through the rotationally broadened 
Ha profile their distances from the stellar rotation axis can 
be measured. These are typically found to be between 2 and 
8 R* and therefore often above the Keplerian co-rotation ra- 
dius of 2.7 R* (Cameron and Robinson 1989a, Cameron and 
Robinson 1989b, Donati and Collier Cameron 1997, Donati 
et al. 1999). These clouds must therefore be supported by 
magnetic loops extending high into the AB Dor corona, and 
are called 'prominences' by analogy to those observed on 
the Sun. Stellar prominences provide us with unique probes 
of the magnetic field topology at the highest levels in the 
stellar corona. 

In Figure[l3] we display the stacked Ha (time-series) 
spectra for each of the four nights of observations. We ap- 
proximate the background in the Ha line by creating a mean 
spectrum that is then subtracted from the entire time-series. 
This leaves only time-varying features such as bright ac- 
tive regions (plage) and rapid absorption transients (promi- 
nences) . When examining Figure[l3] the most striking ob- 
servation is the lack of prominences seen at this epoch. 
Despite the fact that almost 90% of the rotation cycle is 
covered over four nights, only four prominences are found. 
This is substantially less than the average number of 6 to 8 
prominences that are observed at most other epochs (Col- 
lier Cameron et al. 2003). Furthermore, the distribution of 
prominences around the star also appears to be peculiar, 
as almost half of the star (between phases 0.3 to 0.8) does 
not appear to possess any prominences that transit the disk, 
and is covered in bright chromospheric emission from plage 
regions instead. This is unusual for AB Dor, as prominences 
are typically distributed at a much wider range of longi- 
tudes. However, a similar asymmetric distribution has been 
reported on K3 dwarf Speedy Mic (Dunstone et al. 2006a). 
Note that the longitudes that are effectively prominence-free 
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Figure 12. A lower temperature model: assuming an isothermal corona with T=2.2MK and log (EM)[cm" 3 ] = 51.8. Left: The rela- 
tionship between the extent of the X-ray emitting corona and the coronal density, assuming the magnetic field map shown in Figure 1 
and the lower temperature coronal model properties. Right: The dependence of filling factor on n e , also for the lower 2.2 MK coronal 
model. The measured density, logn e [cm — 3 ] = 10.3, suggests a very compact corona with a very similar source surface radius value to 
that estimated in the case of the 10 MK corona, R ss ~ 1.3R* compared to 1.4iJ* for the hotter corona. The filling factor for this model 
would be slightly higher: / = 0.15 compared to 0.11. 



(phases 0.3-0.8) correspond to the most magnetically-active 
longitude region in the surface spot and magnetic field maps 
(compare Figure[l3] with Figure|3J. 

We use the simple technique of Gaussian-fitting (Col- 
lier Cameron 1989a) to track the Ha absorption signatures 
and measure the spatial and radial positions of the corre- 
sponding prominences. Our best-fits are shown in Figure [T3l 
and our measurements of the phases and distances of these 
prominences are displayed in Table [5] Despite prominences 
being detected at similar phases on successive nights their 
heights vary a lot indicating that they may not be exactly 
the same structures from night to night. Instead they are sta- 
ble locations in which prominences tend to reform quickly 
after being ejected. Further details and analysis of these 
Ha observations will be described by Dunstone et al. (in 
prep.). Prominences A, B and C are at distances of between 
2-4 R*, and are typical of prominences found on AB Dor. 
Prominence D is quite different as it is at a high latitude 
and must be relatively close to the stellar surface to tran- 
sit the disk: the distance from the rotation axis indicated 
by its drift rate is 0.91 R* on December 11 and 0.52 R* on 
December 12. This prominence first appears as an absorp- 
tion feature in the Ha line at a velocity of -10 kms -1 as 
opposed to the ~-60 kms" 1 of the other prominences. This 
confirms that the prominence starts to transit the stellar 
disc at a high latitude (~60°). Unfortunately as the phase 
coverage is poor and the end of prominence D's transit is 
not observed our determination of its exact height is uncer- 
tain. It is nevertheless unusual to find such a low lying, high 
latitude, prominence on AB Dor. 

In addition to the prominence absorption features we 
also track two bright plage regions, labelled PL1 and PL2 
in Figure 1131 Both of these, like prominence D, are found 
to have distances from the rotation axis less than 1 _R*. 
If we assume that PL1 and PL2 are indeed active plage 
regions on the AB Dor surface then their latitudes are 44° 
and 49° respectively. It is also interesting how near their 
phases of observation are to those of prominences B and C 



respectively. This is similar to the observation by Jeffries 
(1993) that prominences on Speedy Mic were followed or 
preceded by strong emission transients. 

Soft X-ray emission from the X-ray coronal active re- 
gions may be absorbed by these transiting stellar promi- 
nences, (e.g., as seen in the X-ray Einstein lightcurve of the 
M5V, Proxima Centauri; Haisch et al. 1983). Now that the 
phases at which prominence material is transiting the stellar 
disc have been measured, we can compare these to the X-ray 
lightcurve. We only consider the parts of the rotation cycles 
over which there is an overlap between the X-ray and optical 
observations as prominences can form quickly, within 2P rot . 
Only prominences A and D are known to exist at the same 
time as the Chandra observation. The possibility of soft X- 
ray emission being shadowed by cool prominence material, 
as observed on the Sun, was briefly discussed in Paper I. 
Stellar prominences on rapidly rotating stars have higher 
hydrogen column densities than their solar namesakes. Col- 
lier Cameron et al. (1990) found a column density for AB 
Dor prominences of Nh ~ 10 20 cm -2 , while more recently 
Dunstone et al. (2006b) found column densities in the range 
Nh =1 - 4xl0 19 cm -2 on Speedy Mic. This is an order of 
magnitude more than the Nh ~ 10 18 cm -2 found for solar 
prominences (e.g. Kucera 1998). 

The hydrogen continuum has an upper bound of 911 A, 
below which the absorption cross-section decreases as A -3 . 
The effect should therefore be strongest in the long wave- 
length region of our data (50-130A). This lightcurve is 
shown in Figure[T4]with the prominence positions overlayed. 
This plot focusses on the prominences that were observed 
strictly simultaneously with the X-ray observations. The 
flux- weighted average wavelength of this lightcurve is 83.9A. 
At this wavelength for hydrogen column densities in the 
range 10 19 to 10 20 cm~ 2 we would expect a prominence to 
absorb between 5 and 45% of EUV photons. If the promi- 
nences also contain about 5% He I, then the absorption in- 
creases to between 15-75%. Furthermore if the prominences 
have a fragmented structure similar to solar prominences, 
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then the local density and corresponding absorption could 
be higher. So there is reason to expect that a significant 
reduction in EUV flux would be observed if a prominence 
eclipsed an EUV emitting region of the corona. In contrast, 
approximately 1-2% of the hard X-ray photons from the 1- 
50A lightcurve are expected to be absorbed; this would not 
be detectable due to the intrinsic variability. 

It should be noted that the only data that are observed 
simultaneously at both X-ray and optical wavelengths lie 
between phases 1.5-2.1, i.e. during the December 11 spec- 
tral time-series. We include the X-ray observations of the 
previous AB Dor rotation. Initial investigation using 1 ksec 
bins (Figure la) revealed a possible weak correlation be- 
tween the location of prominence D and a sharp dip in the 
X-ray light curve. To establish the duration of this event 
more accurately we re-binned the light curve into 500 s bins 
to produce Figure [T4l The dip is stronger but only appar- 
ent in one data point; we therefore increased the sampling 
further to produce a lightcurve with 250 s bins and found 
that this point split into two equal fluxes, which lay below 
the other datapoints. So it appears that for 500 s the EUV 
flux decreases considerably and it does so at a phase that 
corresponds closely to that of the central meridian crossing 
of prominence D. The phase difference is A<f> = 0.007 (or 
~ 300 s) , which is the same as the exposure time of each Ha 
spectrum. 

In order to establish the strength of this 500 s dip and 
its significance we need to compare it to a reference flux level 
to mark the unmodulated emission level of the lightcurve. 
We measure this three ways: (a) the local average (average 
of the counts of the four data points immediately to either 
side of the dip), (b) the optimal average flux of the whole 
light curve, and (c) the minimum background level (found 
by taking an iterative 10<r clip of all points above the opti- 
mal average level). The local and optimal average flux level 
are similar: 134 and 139 counts respectively; the magnitude 
of the dip is therefore between 52-54% corresponding to a 
significance of 8— 9<r. The minimum background level is 114 
counts, indicating a dip of 43% and a significance of 6a. Note 
that the significance of the dip discussed here only corre- 
sponds to the error associated with the flux measurement; 
there is also added uncertainty introduced by the intrinsic 
short-term variability observed in the Chandra lightcurves. 
The short wavelength (1-50 A) lightcurve shows no signif- 
icant dip at the corresponding phase (see Figure la). This 
lightcurve shows variations consistent with the intrinsic vari- 
ability of the star; with variability around this phase being 
between 1-2<t (depending on exactly how the mean level of 
emission is defined). This is consistent with the scenario of 
a cool prominence-type structure causing the significant dip 
observed in the EUV lightcurve as the effect of this type 
of cool structure on the 1-50A lightcurve is expected to be 
much smaller (between 1-2% compared to 15-75% in the 
EUV lightcurve). 



7 DISCUSSION 

7.1 Correlation between active longitudes at 
different types of magnetic activity 

There is a significant level of correlation in the positions of 
active longitudes when comparing the surface magnetic field 
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Figure 14. EUV light curve obtained by summing up counts over 
the 50 - 130 A range, focussing on the phases over which there are 
simultaneous optical and X-ray observations. The hashed regions 
show the phases at which different prominences are transitting 
the stellar disc, while the vertical black dashed lines show the 
phase of central meridian crossing. Note the correlation between 
a rapid dip in the lightcurve at phase 2.03 and prominence D. 



maps with spot maps as well as when comparing magnetic 
field maps or spot filling factors with the X-ray lightcurve. 
Some agreement has been observed between the positions of 
active regions as traced in surface spot and magnetic field 
maps of AB Dor at previous epochs, although this is not 
found in similar studies of other systems (e.g., K0V, LQ Hya, 
or the primary component of the RSCVn binary, HR1099) 
(see Donati et al. 2003). The correlation may not be so ap- 
parent in studies of these other systems as they are more 
slowly rotating and it is therefore more difficult to obtain 
the same degree of phase coverage and hence detail in their 
surface maps. Further analysis of detailed surface magnetic 
field maps of stars at different epochs would be necessary to 
establish how common this degree of correlation really is. We 
have carried out a similar analysis of spot and magnetic field 
maps of AB Dor obtained at previous epochs (1998 January, 
1996 December and 1995 December) and do find some cor- 
relation between spot and magnetic field active longitudes 
although at a much weaker level than at 2002 December (the 
epoch analysed in this paper). Hence, it appears that 2002 
December is an unusual epoch for AB Dor in two respects: 
(a) half of the star, between phases ~ 0.1-0.5, appears to be 
covered in stronger magnetic fields, increased spot coverage, 
and stronger chromospheric emission compared to the other 
half; and (b) cool Ha absorbing prominences are not dis- 
tributed across all longitudes as in previous epochs, rather 
they appear to be restricted to the less active phases (see 
FigurerT3J. At this epoch there appears to be a clear link 
between activity at all observed atmospheric levels. 

This is the first dataset for which surface spot, magnetic 
field maps have been acquired contemporaneously with X- 
ray observations and we clearly find significant correlation 
between both surface active longitudes (as traced by the spot 
and magnetic field maps) and the coronal active longitudes 
(as traced by the X-ray lightcurve). This correlation would, 
of course, have been obscured if there had been any large 
flares on AB Dor over the course of the X-ray observation. 
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Figure 13. Raw Ha time series spectra for 2002 December 11, 12, 13 & 14, with phase plotted against velocity. Dashed black lines show 
the stellar v e sin i limits. The grey-scale runs from black at 0.80 times the local continuum level, to white at 1.05 times continuum. The 
dark absorption transients indicate the presence of cool prominences lying above the stellar surface. H« emission traces the locations 
of bright plage regions as they transit the projected stellar disk. By measuring the rate of these transits the heights of the absorbing 
and emitting structures can be deduced. A mean spectrum has been subtracted from the entire time-series. For reference with Table 4 
prominences are labelled A,B,C and D and two of the bright surface regions are labelled PL1 and PL2. 



Table 5. The results of the prominence tracking analysis displayed with phases of meridian crossing and calculated distances from the 
rotation axis (-]£-)■ 



Feature December 11 December 12 December 13 December H 
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Phase 
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A 


0.913 


2.64 ± 0.05 


0.910 


2.55 ± 0.03 


0.896 


3.48 ± 0.05 


0.887 


3.84 ±0.12 


B 






0.817 


3.74 ±0.04 


0.811 


2.55 ±0.03 


0.819 


2.88 ±0.03 


C 










0.247 


2.00 ±0.04 


0.231 


3.31 ±0.05 
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0.039 


0.91 ±0.03 


0.035 


0.52 ±0.19 










Plage 1 


0.809 


0.72 ±0.02 
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0.274 


0.66 ±0.02 











We investigate what level of correspondence we can realis- 
tically expect between surface and X-ray active longitudes 
using coronal models based on extrapolations of the surface 
magnetic field map of AB Dor in Section 4.1 and find that 
the unobserved field in the obscured hemisphere can have a 
significant impact on the shape of the X-ray lightcurve (see 
Figures 7 and 8). Hence, we cannot accurately reproduce 
the exact shape of the observed X-ray lightcurve without 
knowing the type of field in the obscured hemisphere. Our 
results suggest that one would expect significant correlation 
between X-ray emission and optical variability caused by 
dark starspots in active main sequence stars. 

7.2 Sizes of X-ray emitting structures compared 
to other coronal features 

We have found that the size of the X-ray corona must be 
very small, with a height that is not greater than 0.4 R*. 
The main factors that limit the size of the X-ray loops are 
(a) high electron density, n e ~ 2 x 10 10 cm~ 3 , and (b) the 
complex multi-polar surface magnetic field distribution that 
is recovered using ZDI. Sanz-Forcada et al. (2003) estimate 
that the 2 MK corona can extend to 1.3 R*, while the 10 MK 



corona would extend out to 0.24 R*. However, their models 
assume a global dipole, for which the magnetic field strength 
decays more slowly with height, r, B ~ r~ 3 ; our ZDI maps 
of AB Dor clearly show a complex multi-polar field distribu- 
tion; which decays much more quickly with height. As X-ray 
loops can only be supported when the magnetic pressure ex- 
ceeds the electron pressure, we find that lowering the tem- 
perature to 2 MK does not change the height of the source 
surface significantly as a compact corona is necessitated by 
the high n e value measured using density-sensitive diagnos- 
tics in the Chandra X-ray data. Our estimate of the coronal 
extent, between 0.3-0.4 R», is consistent with the 0.3 R* flar- 
ing loop height determined from a large X-ray flare observed 
on AB Dor by Maggio et al. (2000). 



Some n e measurements of the hotter, 10 MK plasma in 
AB Dor and other active cool stars indicate densities that 
are even higher than that measured at 2 MK with 12 < 
logn e [cm -3 ] < 13 values (e.g., Sanz-Forcada et al. 2002, 
Testa et al. 2004). This would imply that a component of 
the corona at 10 MK is extremely compact, with a height, 
H < 0.1 R„. 
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7.3 Complex multi-thermal coronae 

We present tentative evidence for EUV absorption due to 
the presence of a cool prominence transiting near the stellar 
pole. This is detected as a strong dip in the EUV lightcurve 
near the meridian crossing of prominence; no corresponding 
dip is observed in the short wavelength (1-50A) lightcurve. 
The short duration of the eclipse would suggest that both 
the prominence and the EUV emitting region are spatially 
narrow in longitude. Given that the observed dip is approxi- 
mately 50 % of the entire EUV flux coming from AB Dor this 
would suggest that the prominence core is very dense and 
narrow and at least one dominant EUV emitting region is 
highly localised and near the pole. This would then indicate 
that at least half of the EUV flux from AB Dor's corona orig- 
inates near the stellar surface and at high latitudes (near the 
stellar pole) in fewer and more concentrated emitting regions 
than in our X-ray models. This region would contribute to 
the unmodulated component of emission in the X-ray and 
EUV lightcurves. However, as there is considerable intrinsic 
variability due to flaring, and the distance and position of 
the prominence could not be tracked fully, we cannot defini- 
tively measure the sizes of the obscuring and EUV emitting 
structures involved. In the future more simultaneous X-ray 
and optical observations would allow us to better understand 
the physical properties of the prominences and the sizes and 
distribution of EUV emitting regions. 

We conclude that the corona is likely to be made up 
of different types of structures, ranging from extremely hot, 
dense (presumably flaring) plasma very near the surface, to 
a less dense corona at a temperature of a few MK extend- 
ing out to 0.3-0.4 R„, leading finally to an extended corona 
with a height that may reach the Keplerian co-rotation ra- 
dius (~2.7R»). This extended corona could support much 
cooler, low density material that would not be detectable in 
X-rays (e.g., cool prominences with temperatures of between 
8000-9000 K). While the models presented here are simple 
X-ray models, allowing only one coronal temperature and 
density; we have demonstrated that these models are use- 
ful in estimating the global properties of stellar coronae and 
can reproduce key characteristics of the X-ray observations: 
the densities, dominant temperatures and rotational mod- 
ulation in the X-ray lightcurves and spectra. Our models 
show that the bulk of X-ray emission observed from the star 
is unmodulated due to a quite large component that is lo- 
cated in the high latitude regions and at the stellar pole 
(see Figure]!)] for X-ray images). In the future, we will build 
more complicated coronal models, allowing a range of tem- 
peratures and densities based on extrapolations of surface 
magnetic field maps. Future multi- wavelength campaigns, 
co-ordinating X-ray, UV observations with Doppler imaging 
and Zeeman Doppler imaging studies of other active cool 
star systems will enable us to understand how representa- 
tive these models are of cool stars in general. 

7.4 Summary 

• We have obtained contemporaneous X-ray and optical 
data probing magnetic activity at the surface and corona of 
AB Dor at the same time. 

• We find significant correlation between the positions 
of active longitudes in all our activity diagnostics as deter- 



mined from the X-ray lightcurve and the surface spot and 
magnetic field maps. At this epoch, AB Dor appears to pos- 
sess one very large active longitude region that displayed 
enhanced spot coverage, magnetic fields and chromospheric 
emission. 

• We produced detailed three-dimensional models of the 
quiescent X-ray corona of AB Dor by extrapolating surface 
magnetic field maps of the system; using coronal densities 
and temperatures determined from X-ray data of the sys- 
tem. Our models indicate that the X-ray emitting corona 
must be concentrated very close to the surface of AB Dor, 
with a height, H ~ 0.3-0.4R,. 

• The X-ray models suggest that the "missing" flux in 
the obscured hemisphere of the star needs to be of the same 
polarity as that in the observed hemisphere in order to re- 
produce the same level of rotational modulation observed in 
the X-ray lightcurves and spectra. 
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